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Cyclic nucleotide-gated (CNG) channels bind cGMP or cAMP in a cytoplasmic ligand-binding domain (BD), and 
this binding typically increases channel open probability (Po) without inducing desensitization. However, the cat- 
fish CNGA2 (fCNGA2) subtype exhibits bimodal agonism, whereby steady-state Po increases with initial cGMP- 
binding events ("pro" action) up to a maximum of 0.4, but decreases with subsequent cGMP-binding events ("con" 
action) occurring at concentrations >3 mM. We sought to clarify if low pro-action efficacy was either necessary or 
sufficient for con action to operate. To find BD residues responsible for con action or low pro-action efficacy or 
both, we constructed chimeric CNG channels: subregions of the fCNGA2 BD were substituted with corresponding 
sequence from the rat CNGA4 BD, which does not support con action. Constructs were expressed in frog oocytes 
and tested by patch clamp of cell-free membranes. For nearly all BD elements, we found at least one construct 
where replacing that element preserved robust con action, with a ratio of steady-state conductances, g(io mM cgmp)/ 
g(3 mM cGMP) < 0.75. When all of the BD sequence C terminal of strand (36 was replaced, g(io mM cGMP)/g(3 mM cgmp) was 
^ increased to 0.95 ± 0.05 (n = 7). However, this apparent attenuation of con action could be explained by an in- 

Q crease in the efficacy of pro action for all agonists, controlled by a conserved "phosphate-binding cassette" motif 

~Q that contacts ligand; this produces high Pq values that are less sensitive to shifts in gating equilibrium. In contrast, 

•Jj; substituting a single valine in the N-terminal helix aA abolished con action (g(3o mM cGmp)/ g(3 mM cGmp) increased to 

1.26 ± 0.24; n = 7) without large increases in pro-action efficacy. Our work dissociates the two functional features of 
n low pro-action efficacy and con action, and moreover identifies a separate structural determinant for each. 

0) INTRODUCTION 
c 

^ The activation of signaling receptors by agonists is ex- 

. , plained mechanistically by classical allosteric coupling 

O theory (Monod et al., 1965): the agonist's interactions 

"fC with its binding site are more stable when the receptor 
E is in an activated conformation compared with in a 
Q deactivated conformation. Negative or inverse agonists 

~~> are ligands that allosterically reduce receptor activity; 
Q> their receptor interactions are "counterproductive," being 

I— more stable for the deactivated conformation than for the 
activated conformation (Costa and Herz, 1989). Efforts 
at efficient rational drug design would benefit greatly 
from a mechanistic understanding of the counterpro- 
ductive interactions that make a ligand specifically into 
a negative agonist instead of a conventional (positive) 
agonist (Bond and Ijzerman, 2006) . Historically, numer- 
ous negative agonists are known for metabotropic (G 
protein-coupled) receptors (Kenakin, 2004) , but there 
are rather fewer known instances of negative agonism in 
ionotropic receptors (ligand-gated ion channels) ; exam- 
ples include the NMDA receptor (Mony et al., 2009) and 
members of the Cys-loop superfamily like the GABAa 
and nicotinic acetylcholine receptors (Olsen et al.. 
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2004; Arias et al., 2006) . (Note that channel pore block- 
ers are widely known, but those work by competitive in- 
hibition rather than by an allosteric mechanism at a site 
away from the channel pore.) 

Recently, a new form of negative agonism has been 
described for cGMP acting on CNG channels. CNG chan- 
nels are nonselective cation channels activated by direct 
binding of cGMP and cAMP; the channels consist of 
homotetramers or heterotetramers of homologous sub- 
units. Each subunit has a conserved architecture of six 
transmembrane helices and a cytoplasmic C-terminal 
region that contains the ligand-binding domain (BD). 
Cells such as photoreceptors and olfactory sensory neu- 
rons express various types of CNG subunits (i.e., para- 
logues) in different combinations (Bradley et al., 2001; 
Kaupp and Seifert, 2002; Matulef and Zagotta, 2003). 
Negative agonism is observed in a particular subtype of 
CNG channels — the CNGA2 subtype from catfish olfac- 
tory neurons (fCNGA2) (Goulding et al., 1992) — and in 
recombinant CNG channels derived from this subtype 
(Young et al., 2001; Chan and Young, 2009). Notably, 
the dose-response relation of activity versus cGMP con- 
centration shows a rising phase for low concentrations and 
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a falling phase for higher concentrations. This contrasts 
with the normal (monotonically rising) dose-response 
curve of the positive agonist cAMP. Single-channel re- 
cordings have shown (Young et al., 2001) that the fall- 
ing phase of the cGMP dose-response reflects a decrease 
in steady-state open probability (Fq) induced by a new 
cGMP molecule binding to a partially liganded channel 
(see reaction scheme in Fig. 1 ) . That is, it does not involve 
a reduction of open-channel unitary conductance as 
typically found with agonist pore block (Ogden and 
Colquhoun, 1985; Karashima et al., 2007), or a sponta- 
neous gate-closing conformational change as found 
with traditional desensitization in ligand-gated channels 
(Cachelin and Colquhoun, 1989; Robert and Howe, 2003) 
or inactivation in voltage-gated channels (Hoshi et al., 
1991). We have termed this unusual phenomenon bi- 
modal agonism because cGMP exhibits two opposite 
modes of action: "pro action" that enhances receptor ac- 
tivity by a positive agonism mechanism, and "con action" 
that suppresses activity by a negative agonism mechanism. 
Our distinct terminology "con action" is used to denote 
negative agonism in this special context where more than 
one agonism mode is possible. 

Besides being specific for cGMP and not cAMP, several 
additional features of con action argue against mecha- 
nisms that involve nonspecific interference with channel 
function. A neutral agonism mechanism, where cGMP is 
displaced from its pro-action site by a competitor mole- 
cule that binds open and closed channels equally, is 
ruled out because the agonist responsible for con action 
is cGMP itself. BD substitution studies showed that al- 
though the fCNGA2 BD imparted bimodal cGMP agonism 
and normal cAMP agonism, the BD from the rat CNGA4 
type (rCNGA4; Bradley et al., 1994; Liman and Buck, 
1994) imparted solely normal cGMP and cAMP agonism. 
Thus, particular sequences within the fCNGA2 BD are 
required for con action (Young et al., 2001; Chan and 
Young, 2009) , as would be expected for a modification 
of gating under BD-mediated control. Negative agonism 
acting through the BD is previously unknown in CNG 
channels because the known "antagonists" such as phos- 
phorothioate analogues (Kramer and Tibbs, 1996) are 
merely weaMy effective ("partial") positive agonists. It may 
be noted that the cGMP affinity of the binding site re- 
sponsible for con action seems to be weak (millimolar) 
in constructs containing the fCNGA2 BD, but it is in fact 
comparable to the cAMP affinity of canonical pro-action 
sites in the same channels. Moreover, the operation of 
con action in heteromeric CNG channels (with fCNGA2 
sequence in only two of four BDs) is conditional on sub- 
unit arrangement, consistent with the notion that the 
mechanism requires a specific geometry in the ligand- 
channel complex (Chan and Young, 2009). 

Bimodal agonism illustrates how it is misleading to say 
a ligand "is" a negative agonist; strictly, we can say only 
that the ligand "acts as" a negative agonist. How a ligand 



acts will depend on the receptor structure, and previ- 
ously bound Ugands must be included in describing that 
receptor structure. In the thermodynamic cycle for pro 
action at the left side of the Fig. 1 reaction scheme, 
cGMP binding preferentially stabilizes the open chan- 
nel only when the number of bound ligands is less than 
some critical number, m. When some partially liganded 
channel state is attained (m ligands in Fig. 1 ) , the cou- 
pling relationship between binding and gating becomes 
reversed. Subsequent binding of the same agonist spe- 
cies now preferentially destabilizes the open state (con 
action) with counterproductive interactions. 

Con action in CNG channels represents a special ex- 
ample of negative agonism that manifests "alongside" pos- 
itive agonism, caused by the same cGMP ligand species 
and relying on the same BD region. Thus, it becomes 
pressing to clarify the causal relationship between them. 
For instance, the previous BD substitution studies found 
that the rCNGA4 BD not only imparted normal cGMP 
agonism but also imparted a much higher efficacy of pro 
action by both cGMP and cAMP compared with the 
fCNGA2 BD. Might the particular structural deficien- 
cies responsible for poor pro-action efficacy in fCNGA2 
be necessary for operation of the con-action mechanism? 
If true, this would imply that any future compounds se- 
lected for effective con action in fCNGA2 will fail to in- 
duce con action in any of the numerous CNG channels 
whose pro-action mechanism is highly efficient. Con- 
versely, might the poor pro-action efficacy be in itself 
sufficient to bring a con-action mechanism into opera- 
tion? If true, this would imply that cGMP-negative ago- 
nism should already operate in every example of a low 
efficacy CNG channel, perhaps going undiscovered be- 
cause higher cGMP concentrations were not previously 
tested. Here, we identify counterexamples against both 
the necessity hypothesis and the sufficiency hypothesis. 
That is, we show that con-action and low pro-action effi- 
cacy, the two distinctive functional features of the fCNGA2 
BD, derive from independently operating structural de- 
terminants in the BD and are thus not causally related. 
Moreover, we isolate a single residue in the fCNGA2 BD 
with a critical role in con action, opening the door for 
future elucidation of a specific mechanism for negative 
agonism in CNG channels. 

MATERIALS AND METHODS 

Recombinant DNA 

Previous reports (Tibbs et al., 1997; Young et al., 2001) have de- 
scribed the design of the "X-chimera" series of BD substitution 
constructs, with invariant non-BD sequence derived from bovine 
CNGAl (bCNGAl; Kaupp et al., 1989) and fCNGA2 (Goulding 
et al., 1992). New X-chimera constructs were generated with pre- 
viously reported PCR/restriction/ligation techniques (Young and 
Ki-Qugliak, 2004), incorporating BD sequences from fCNGA2, 
residues L455-A583, and from rCNGA4 (Bradley et al., 1994; 
Liman and Buck, 1994), residues L356-A484. Some constructs 
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were previously reported with different names in Young et al. 
(2001): X-fA2 is synonymous with X-a, X-rA4 with X-p, Construct 
1 with X-o[r/Pc, and Construct 6 with X-Pr/olc; in addition, Con- 
struct 2 is synonymous with X-bimP reported previously (Chan 
and Young, 2009). There are six constructs new in this study 
(not including point mutants), with the following concatena- 
tions: Construct 3, fCNGA2 (L455-C511) + rCNGA4 (;F413-A484); 
Construct 4, fCNGA2 (L455-K488) + rCNGA4 (E390-A484); Con- 
struct 5, fCNGA2 (L455-K488) + rCNGA4 (E390-S426) + fCNGA2 
(G526-A554) + rCNGA4 (E456-A484); Construct 7, rCNGA4 
(L356-R389) + fCNGA2 (E489-A583); Construct 8, rCNGA4 
(L356-G372) + fCNGA2 (D472-A583); Construct 9, fCNGA2 
(L455-G471) + rCNGA4 (E373-R389) + fCNGA2 (E489-A583). 
Sequence regions subjected to PCR were confirmed by dide- 
oxy sequencing. 

Electrophysiology 

Oocytes were harvested from mature female Xenopus laevis frogs 
(procedures approved by the Canadian Council on Animal Care) , 
injected with in vitro-transcribed RNA, and assayed in the excised 
inside-out patch clamp configuration, with all procedures done as 
previously reported (Young et al., 2001; Young and Krougliak, 
2004; Chan and Young, 2009). Recording conditions are briefly 
summarized: Pipette and bath solutions both contained (in mM): 
67 KCl, 30 NaCl, 10 HEPES, 10 EGTA, and 1 EDTA, pH 7.2 with 
KOH. Na salts of cAMP or cGMP were included by iso-osmolar re- 
placement of NaCl. For Ni^* potentiation experiments (Young et al., 
2001), the bath solution was modified to include 10 pM NiClj 
(Sigma-Aldrich) , omit EGTA and EDTA, and increase KCl to 89 mM. 
Current measurements were acquired at room temperature 
(19-23°C) with either an Axopatch 200B or GeneClamp500 am- 
plifier, filtered at 10 kHz using the amplifier's low-pass Bessel fil- 
ter, and recorded at 0.1-1 kHz by a Digidata 1322 and pClamp 9.0 
software (Axon Instruments) . 

Patches were held at —40 mV, and steady-state currents in the 
presence of cyclic nucleotide were corrected by subtraction of 
leak currents in the absence of cyclic nucleotide. Spontaneous ac- 
tivity changes (run-up and run-down) of CNG channels in excised 
patches were typically complete after 10-20 min, somewhat lon- 
ger than previously noted for bCNGAl (Molokanova et al., 1997) . 
With the exception of specific experiments examining this run-up 
(Fig. 3 B), all reported measurements were collected after run-up 
was completed, as shown by a difference of <10% among multiple 
conductance measurements in 3 mM cGMP. Steady-state conduc- 
tance in 3 mM cGMP after completion of run-up was used in each 
patch to normalize conductances in other conditions, including 
those before the completion of run-up (Fig. 3 B) . Unless otherwise 
noted, means are reported ± SD with n the sample size, and con- 
ductance ratios were compared with unity using t test. 

For selected constructs, the Po for the 3-mM cGMP normaliza- 
tion condition (after run-up) was estimated as the conductance ratio 
g(3 mM ccuT) / g(m mM cAMP+Ni), where the denominator is the conduc- 
tance in the presence of 30 mM cAMP and 10 pM Ni^* ion 
(Varnum et al., 1995; Young et al., 2001). Ratio values recorded for 
constructs in this study were as follows: Construct 2, 0.180 ± 0.050 
( n = 3) ; Construct 3,0.73±0.15(re=3); Construct 8, 0.337 ± 0.092 
(n = 4); X-fA2 V457E, 0.359 ± 0.079 {n = 4). This method as- 
sumes Po » 1 in the nickel-potentiated fully liganded channel and 
that unitary open-channel conductance is unaffected by either 
cyclic nucleotide or Ni^*. These two assumptions were vali- 
dated for previously studied X-chimeras (Young et al., 2001); all 
X-chimera sequences are identical in the transmembrane domain 
(including the pore region determining unitary conductance; 
Goulding et al., 1993) and the C-linker (including the Ni^*-binding 
histidine residue; Gordon and Zagotta, 1995). Dose-response 
relations in Figs. 4 and 5 were fitted with the Hill equation, P^ = 

Pmax/(l+Kl/2/[cNMP])*. 



Comparative models 

Models were built for the X-fA^2 C-linker and BD (211 residues, 
consisting of bCNGAl M401-G484 followed by fCNGA2 L455- 
V581) and for several constructs and mutant forms studied. Each 
target sequence was aligned to a single template of D443-L643 
from mouse HCN2 (Protein Data Bank accession no. 1Q50; 
Zagotta et al., 2003), which was the closest relative with struc- 
tural data available; there was 27% identity with three gaps intro- 
duced in the template sequence. Models were built using the 
"Fit Raw Sequence" function in DeepView 4.0.1 (Guex and 
Peitsch, 1997) and the Project mode of SWISS-MODEL (version 
8.05) (Arnold etal., 2006). Root-mean-squared distance calcula- 
tions and detection of atoms within 5 A of particular side chains 
were performed using DeepView; hydrogen bonds were detected 
using VADAR (version 1.6) (Willard etal., 2003). For model images 
(Fig. 4 C), PYMOL (DeLano Scientific) was used for the ribbon 
diagram and GRASP2 for the electrostatic potential (Petrey and 
Ho'nig, 2003). 

RESULTS 

Rationale of BD substitution 

The example macroscopic current trace in Fig. 1 illustrates 
pro and con action in homomeric bimodal CNG chan- 
nels formed by a recombinant chimeric construct, X-fA2, 
whose composition and functional properties have been 
described in detail (Young et al., 2001). X-f\2 contains 
the fCNGA2 BD sequence and reconstitutes the gating 
and permeation characteristics of intact fCNGA2, just 
with improved levels of homologous expression resulting 
from selected non-BD sequence substitutions. Notably, 
the dose-response curves of X-£\2 for both cGMP and 
cAMP match those of intact fCNGA2, including matching 
maximal ?„ values (Young et al., 2001). A subsequent 
series of BD substitution constructs (X-chimeras) were 
bcised on substituting BD sequence in the original X-fA2 
parent (Young et al., 2001; Young and Krougliak, 2004; 
Chan and Young, 2009) . This enables comparisons of the 
activation properties of various BD sequences, set within 
one context of invariant non-BD sequence (Fig. 2 A) . 
Each X-chimera subunit can be expressed alone in 
Xenopus oocytes to form functional homomeric CNG 
channels, which are then assayed as in Fig. 1 (inside-out 
excised patch clamp with agonist solutions perfused on 
the cytoplasmic side). 

Perfusion of 3 mM cGMP on X-fA2 homomers elicits 
a nondesensitizing channel current (pro action; Fig. 1 
trace) , and throughout this study we use the steady-state 
conductance elicited by 3 mM cGMP, g^^ „,m cgmp)> as a 
normalization reference for steady-state conductances 
in other agonist conditions in the same patch (see Fig. 1 
caption and Young et al., 2001 for interpretation of 
"spiked" current traces showing con action with cGMP 
dosages >3 mM) . The BD incorporated in each X-chimera 
tested can then be classed as either bimodal or normal, 
according to whether the cGMP dose-response curve of 
the X-chimera exhibits the peculiar negative slope asso- 
ciated with con action. 
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Within the 130-residue BD sequence of CNG chan- 
nels, there are 36 residues that are not conserved between 
fCNGA2 and rCNGA4 (Fig. 2, A, alignment, and B, 
black ticks) . In this study, we substituted subsets of these 
36 nonconserved residues with their rCNGA4 counter- 
parts and tested the effects on cGMP gating. The BD is ho- 
mologous to those of other cyclic nucleotide-dependent 
proteins (Weber et al., 1982; Kaupp et al., 1989), with an 
N-terminal helix (ctA), a (3-roll fold (strands (31-(38), and 
two C-terminal helices (aB-C); an additional "PB" helix 
and loop interrupt the (3 roll. The PB elements form a 
conserved "cassette" that participates in ligand binding 
in all known cyclic nucleotide-dependent proteins using 
the BD (Diller et al., 2001). The 36 nonconserved resi- 
dues are spread throughout the BD. Any of the fCNGA2 
BD residues that we substituted might be "critical" for 
con action, meaning that its substitution would abolish 
con action; there might be several such critical residues. 
Or, substitution of one or more fCNGA2 residues might 
not abolish con action but might instead modify quanti- 
tative features of the falling phase in the cGMP dose- 
response curve. Such residues would be properly described 
as "contributing" to con action but not critical for it. 
In addition, substitutions might change pro-action effi- 
cacy; this would modify the quantitative features of the 
rising phase of the dose-response curve. The key ques- 
tions are whether it might be possible to increase pro- 
action efficacy without abolishing con action (showing 
that low pro action is not necessary for con action) , and 
whether it might be possible to abolish con action with- 
out increasing pro-action efficacy (showing that low pro 
action is not sufficient for con action) . 

Extensive BD sequence substitutions preserve con action 
We had previously constructed Constructs 1 and 2 with 
portions of the C-terminal region of the fCNGA2 BD 
replaced with rCNGA4 residues (Fig. 2 B) ; both of these 
constructs were bimodal (Young et al., 2001; Chan 
and Young, 2009) . We extended the X-chimera series 
with Constructs 3 and 4, replacing larger amounts of 
C-terminal sequence from the fCNGA2 BD, and these 
progressive replacements increaised the value of the ratio 

^■(10 mM cGMP)/ giS mM cGMP)- Indeed, the g^w mM cGMP)/ g{3 mM 

cGMP) value for Construct 3 is very close to unity (although 
still significantly less, P < 0.02), and the value for Con- 
struct 4 is not significantly different than unity (P > 0.3) . 
These results might suggest a critical role for some of 
the 10 fCNGA2-specific residues from (35 through the PB 
cassette. However, surprisingly, we were able to con- 
struct a counterexample to this hypothesis: starting 
from a Chimera 4 background, we added back a set of 
seven fCNGA2-specific residues in another region, |37 
through the B helix, and in doing so we "rescued" con 
action. The resulting Construct 5 has a ratio of g^io mM 
cGMP)/g"(3 mM cGMP) = 0.72 ± 0.12 (w = 7), significandy below 
unity (P < 0.001). 



Considering the two bimodal Constructs 2 and 5 to- 
gether enables us to exclude a critical role for a large 
number of residues that are unique to the fCNGA2 
BD. Specifically, each unconserved residue C-terminal 
of strand (33 has been replaced in either Construct 2 
or Construct 5, or both, without abolishing con ac- 
tion. On the contrary, con action is robustly main- 
tained with g-(io mM cGMP)/ g(3 mM cGMP) < 0.75. Therefore, 
if there exists a critical residue, it must be located at 
the amino end of the BD, within the region spanning 
helix aA through strand (33 (to be discussed further in 
a later section) . 



3mM 10 mM 

cGMP cGMP 




cGMP cGMP cGMP 

^ Con action 

Pro action 



Figure 1 . Simple allosteric reaction scheme for bimodal ago- 
nism. Reaction scheme (Chan and Young, 2009) is illustrated by a 
cell-free inside-out macroscopic patch recording of bimodal CNG 
channels with cGMP washed into and out from the bath (solid 
and dashed black bars, respectively). Applying cGMP elicits in- 
ward channel current with no desensitization up to 3 mM cGMP, 
but with 10 mM cGMP, the steady-state macroscopic conductance 
(white arrow) is smaller than for 3 mM. During wash-in and wash- 
out of the 10-mM cGMP solution, there are brief intervals with 
bath cGMP concentration near 3 mM giving rise to "spikes" of 
maximal conductance. Details of reaction scheme notation: 
Closed- (C) and open- (O) channel states are equilibrating rap- 
idly compared with agonist wash-in and wash-out. Sizes of black 
reaction arrows indicate the favored direction of equilibrium. 
Po increases with cGMP binding of up to m molecules of cGMP 
(pro action). Binding of the (m-Fl)* cGMP molecule decreases 
the steady-state Po (con action) to below the value obtained with 
m ligands. Experimental details of current trace: The previously 
characterized channel, X-fA2, incorporates the BD from fCNGA2 
channel in a chimera with sequence from other CNG channel 
types (see Young et al., 2001, and Fig. 2). X-fA2 was expressed 
as homomers in Xenopus oocytes; patches were held at —40 mV. 
White arrows mark times where steady-state conductance, g, was 
measured with cGMP concentrations fixed at their nominal con- 
centrations; for this example, g-(io min cGMP)/g'(3 min cgmp) = 0.63. 
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Substitution of the fCNGA2 PB cassette increases 
pro-action efficacy 

We selected Construct 3 for a more detailed compari- 
son with the previously characterized Construct 2, because 
these two constructs differ only in five residues within 
the key PB cassette region. In multiple cyclic nucleotide- 
activated proteins, movement of the PB cassette in re- 
sponse to ligand binding is believed to be essential for 
propagating conformational change to other parts of the 
BD (Canaves and Taylor, 2002; Rehmann et al., 2003; 
Clayton et al., 2004; Kim et al., 2005; Mazhab-Jafari 
et al., 2007). Thus, although we knew that the five 



fCNGA2-specific residues in the PB cassette were not 
critical for con action, we anticipated that they might be 
important determinants of the low pro-action efficacy 
ofX-£\2. 

We converted macroscopic conductances into Po values 
by determining the maximum conductance in each patch 
(Po = 1) elicited by 30 mM cAMP and 10 pM Ni^^ (see 
Materials and methods) . The cGMP dose-response curves 
(Fig. 3 A, down triangles vs. gray solid curve) show that 
both Constructs 2 and 3 are bimodal with a con phase at 
concentrations >3 mM. However, the maiximal Po at 3 mM 
cGMP for Construct 3 (0.73 ± 0.15; n = 3) is far higher 
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Figure 2. Extensive BD sequence substitutions preserve con action. (A) Topology of the CNG channel subunit, highlighting the BD 
region substituted in the X-chimera series of constructs. Alignment compares the substituted BD sequences (bounded by vertical dotted 
lines) originating from rCNGA4 (in X-rA4) and fCNGA2 (in X-fA2). Invariant non-BD sequences (Young et al., 2001) include bCNGAl 
sequence in the G-linker and the extreme C-terminal region, respectively, before and after the BD. Dots in the sequence alignment 
indicate fCNGA2 residues conserved with rCNGA4; secondary structure elements are marked (oi for helices, (3 for strands) as predicted 
by comparative modeling (Fig. 4) . (B) X-chimera testing substitutions in the C-terminal portion of the BD. Bar represents BD sequence; 
below the bars, secondaiy structure elements are marked (letters for helices, numbers for strands, and PB for the PB cassette) . Gray back- 
ground color in bar represents amino acids conserved between fCNGA2 and rCNGA4; black ticks mark unconserved amino acids found 
in fCNGA2 but not in rCNGA4. Mean g^w mu cGMP)/g(3 mM cgmp) values (±SD, number of patches in parentheses) include those previously 
reported for Construct 1 (Young et al., 2001) and for Construct 2 and X-rA4 (Chan and Young, 2009). (C) Examples of current traces 
for selected constructs during 10-mM cGMP pulses, with g-(io mm cGMF)/g{3 mM cgmp) as marked (based on applications of 3 mM cGMP; not 
depicted). Horizontal bar, 10 s; vertical bar, 1 nA. 
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than for Construct 2 (0.18 ± 0.05; n = 3) , so that the entire 
cGMP dose-response curves for the two channels are 
dramatically separated. Thus, introducing the rCNGA4 
PB cassette sequence into the background of Construct 2 
resulted in higher Fq values for cGMP agonism. How- 
ever, with a bimodal dose-response curve, cGMP pro 
action cannot be reliably quantified separately from con 
action; even at low cGMP concentrations, some binding 
events may be operating with con action at the same 
time as other binding events are operating with pro ac- 
tion. Thus, we entertained an alternative explanation 
where the PB cassette substitution in Construct 3 merely 
attenuated cGMP con action without affecting pro action. 

To clarify the effect of the PB cassette on pro action, 
we evaluated efficacy for cAMP, which appears never to 
participate in the con-action mechanism. The Construct 
2 background channel had mean Pq values ranging firom 
0.36 at 3 mM cAMP to 0.62 at 30 mM cAMP; at compa- 
rable concentrations, the Po for Construct 3 was higher, 
with means ranging from 0.82 at 3 mM cAMP to 0.94 at 
30 mM cAMP (Fig. 3 A, up-triangles vs. gray daished curve) . 
The trajectory of PB cassette movement in the normal 
(pro) action of cyclic nucleotide-activated proteins is 
proposed to be essentially the same whether the bind- 
ing agonist is cGMP or cAMP (Canaves and Taylor, 2002; 
Flynn et al., 2007) . Thus, our results are consistent with 
PB cassette replacement in Construct 3 causing a gen- 
eral increase in efficacy for both cAMP and cGMP. The 
highly efficacious function of the rCNGA4 BD that was 
previously reported (Young et al., 2001) is also consis- 
tent with associating high efficacy with the rCNGA4 PB 
cassette sequence. 

Necessity hypothesis: construct 3 exhibits efficient con 
action that is masked by its high pro-action efficacy 

In previous examples of bimodal CNG channels such 
as fCNGA2, X-fA2, and Constructs 1 and 2, the P^ elic- 
ited by cGMP was relatively low, limited to 0.4 at most. 
This might suggest that some defect in the structural 
apparatus for pro action might be necessary (a pre- 
requisite) for con action. This necessity hypothesis 
could be disproven if a counterexample were found, 
namely a channel with both efficient con action and 
efficient pro action. In this light. Construct 3 provides 
the first example of a bimodal CNG channel whose 
cGMP pro action is efficient enough to produce Po > 0.7. 
We note that Construct 3 also had a higher value for 
^■(10 mM cGMP)/g'(3 mM cGMP) Compared with the other 
bimodal channels, which might suggest that its con- 
action mechanism was not as efficient. However, as 
we explain below, it is not true that an increase in g-(io mM 
cGMP)/g'(3 mM cGMP) always indicates attenuation of con 
action in mechanistic or energetic terms, and the dose- 
response curves of Construct 3 match what would be 
predicted if the con action was unaffected by PB cas- 
sette replacement. 



The Po parameter traditionally plotted in dose- 
response curves does not bear a simple relation with the 
underlying energetics of gate opening; the more perti- 
nent parameter is the equilibrium constant for the pore- 
gating reaction ([open] /[closed], in a simple model of 
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Figure 3. PB cassette substitution masks con action by raising gat- 
ing efficacy: refutation of necessity hypothesis. (A) Dose-response 
relations of Construct 3 activated with cGMP (down-triangles) or 
cAMP (up-triangles). Points show means and error bars show SD 
(n > 3 for each point); points are joined by straight lines without 
fitting. To plot Po values, conductances were normalized to the 
steady-state conductance measured at 3 mM cGMP (solid down- 
triangle) in the same patch. Then, Po at 3 mM cGMP was fixed as 
0.73 based on Ni^* potentiation (see Materials and methods). For 
comparison, gray curves show dose-response relations for Con- 
struct 2 activated by cGMP (solid) or cAMP (dashed); these were 
previously reported on a relative scale without conversion to Po 
values (Chan and Young, 2009), and new Ni^* potentiation exper- 
iments in this study fix the Po at 3 mM cGMP as 0.18. All measures 
for both constructs in this plot were collected after completion 
of spontaneous run-up. (B) Example macroscopic current trace 
from one patch of Construct 3 before and after completion of 
run-up, tracking progressive increase in pro-action efficacy and 
masking of con action. Po scale at right is based on fixing Po = 0.73 
for 3 mM cGMP after run-up. Arrows and asterisks, respectively, 
mark steady-state and spike currents as in Fig. 1. For this example, 
Psteady/Pspiie = 0.80 before run-up and 0.95 after run-up. 
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a single open and closed state). As the numerical value 
of Po increases toward unity, the Po parameter becomes 
less sensitive to changes in the gating equilibrium constant 
(Varnum and Zagotta, 1996). Thus, if a structural modi- 
fication (for example, PB cassette replacement) were 
to increase pro-action efficacy without affecting other 
gating features, this would as a consequence also flatten 
the Pq versus concentration curve, so that significant 
concentration-dependent changes in gating energetics 
caused by either pro or con action would be "masked." 

The gating equilibrium constant cannot be evaluated 
with precision at a subsaturating agonist concentration 
because the individual channel molecules in the popu- 
lation will have different binding site occupancies and 
hence different gating equilibrium constants with differ- 
ent corresponding open probabilities. As a surrogate pa- 
rameter, we can calculate Z = Po/(l— Po) from the 
average Po at a given subsaturating cGMP concentra- 
tion; Z represents an approximate central value, around 
which the individual equilibrium constants in the chan- 
nel population will be clustered. This analysis shows that as 
cGMP concentration increases from 3 to 10 mM, Z de- 
creases by a similar factor for both Constructs 2 and 3. 
For Construct 2, Z changes by a factor of 0.7, decreasing 
from 0.23 at 3 mM cGMP (Po = 0.18) to 0.15 at 10 mM 
cGMP (Po = 0.13); for Construct 3, Z changes by a factor 
of 0.8, decreasing from 2.6 at 3 mM cGMP (Po = 0.73) to 
2.2 at 10 mM cGMP (Po = 0.69) . Although we emphasize 
that Z is not numerically equivalent to a microscopic equi- 
librium constant, the similarity of Z shifts in Constructs 

2 and 3 is consistent with a scenario where the energet- 
ics of con action are similar in these two channels. 

Further evidence for the "masking" of con action in 
Construct 3 is observable during spontaneous "run-up" 
of channel activity in the minutes immediately after 
patch excision. Run-up in X-chimeras, as in many CNG 
channels (Molokanova et aJ., 1997) , typically causes a po- 
tentiation of cyclic nucleotide activation, including an 
increase in efficacy (unpublished data). For patches of 
Construct 3, we noticed that waish-in and waish-out current 
spikes indicative of bimodal agonism were more promi- 
nent in pulses of 10 mM cGMP given early after patch 
excision compared with later pulses (Fig. 3 B) . To quan- 
tify con action in individual pulses of 10-mM cGMP per- 
fusion, we compared the Po value at steady state (Psteady) 
and at the peak of the wash-in spike (Pspike) , which cor- 
responds to an interval where cGMP is momentarily near 

3 mM (Fig. 1); this avoids reference to the steady-state 
g{3 mM cGMP) Standard, which would not be stable during 
run-up. In four patches of Construct 3, pulses of 10 mM 
cGMP were applied less than 5 min after excision (mean 
Psteady = 0.53 ± 0.02) as wcll as after completion of run-up 
(Psteady = 0.69 ± 0.04) . Cou actiou in these four patches was 
more robustly exhibited before run-up (Psteady/Pspike = 
0.84 ± 0.04) compared with after (Psteady/Pspike = 0.97 ± 
0.03; significandy different than before by paired t test. 



P < 0.008). This is consistent with the masking of con 
action arising progressively as Pq increaises during nm-up. 

In summary, the difference in dose-responses of Con- 
structs 2 and 3 can be explained as a potentiation of pro 
action that is general (i.e., applying to both cAMP and 
cGMP) ; the con phase in the cGMP dose-response of 
Construct 3 appears flattened not because of attenua- 
tion of cGMP con-action energetics but because of high 
Po masking. Construct 3 possesses both efficient con ac- 
tion and efficient pro action, which suggests strongly that 
the poor pro-action efficacy exemplified by X-fA2 is not 
a prerequisite mechanistic condition for the operation 
of con action. That is, low pro-action efficacy is not nec- 
essary for con action. Construct 3 does not exclude the 
possibility that a disabled pro-action mechanism could 
be sufficient in itself to induce operation of con action. 

Sufficiency hypothesis: abolishing con action without high 
pro-action efficacy 

Masking with high pro-action efficacy represents one way 
that the manifestation of con action can be suppressed, 
but it does not represent a way to disable the operation 
of the con-action mechanism. We found above (Fig. 2) 
that many fCNGA2 residues from p5 through aC were 
not critical for con action. Therefore, we performed an- 
other mutagenesis search with progressive replacements 
in the N-terminal region of the BD. We knew from pre- 
vious work (Young et al., 2001) that con action could in- 
deed be abolished by replacing very extensive N-terminal 
sequence (Fig. 4 A, Construct 6). The new Constructs 7 
and 8 contain BD sequences >94% identical to the 
fCNGA2 BD, yet like Construct 6, they both lack con ac- 
tion (Fig. 4 A) ; Wcish-in and wash-out current spikes were 
never observed either before or after run-up, and steady- 
state activity increased with cGMP concentration up to 
30 mM (Fig. 4 B, solid down-triangles). Thus, substitu- 
tion of as few as three residues (in helix aA and strand 
(31 of Construct 8) abolished observable evidence for 
con action in the range of 3-30 mM cGMP (we phrase 
this statement conservatively because we could not test 
cGMP concentrations above 30 mM because of viscosity 
and solubility problems, and it is in any case impossible 
to prove a negative assertion for all possible conditions) . 
Notably, bimodal agonism was observed in another Con- 
struct 9 where strands |32 and strands (33 were replaced, 
showing that those sequence elements are not critical 
for con action. 

Although the cGMP response relation of Construct 8 
is not suppressed by a con phase, there is no evidence for 
high Po that could mask con action as in Construct 3. 
The Po for Construct 8 at 3 mM cGMP was only 0.37 ± 
0.09 (n = 4) and did not exceed 0.5, even at 30 mM 
cGMP For both cGMP and cAMP, the Construct 8 dose- 
response curves are very similar to those of X-fA2 (Fig. 4 B, 
black points vs. gray curves) , with the exception of cGMP 
>3 mM. This shows that the c(A-(31 replacement did not 
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observably potentiate pro action for either agonist, a quite 
different effect than that of PB cassette replacement. 

The lack of con action and the low pro-action efficacy 
in Construct 8 refutes the hypothesis that a low pro- 
action efficacy is sufficient to produce a con-action mech- 
anism. Together with Construct 3, our findings show that 
low pro-action efficacy is neither necessary nor sufficient 
to cause the con action found in X-fA2. This strongly 
suggests that there is no causal relation between con 
action and low pro-action efficacy as they are found in 
X-fA2 and intact fCNGA2. 

V457 in helix aA rather than R464 in strand pi is critical 
for bimodal agonism 

Might the critical aA-(31 region contribute to a distinctive 
cGMP-binding site responsible for con action? There is 
no available high resolution structural information on 
the BD of CNG channels, so we assessed the approximate 



location of the key aA-(31 region in the fCNGA2 BD 
using comparative modeling (see Materials and methods) 
based on the crystal structure of a C-terminal fragment 
from the related HCN2 channel (Zagotta et al., 2003). 
Models were built for the C-linker and BD of X-fA2 and 
Constructs 7 and 8 (Fig. 4 C) ; the independently built 
models were closely superimposable with each other 
with a small root-mean-squared distance (0.2 A) for all 
atoms excluding unconserved residues. In other words, 
introducing up to seven substitutions to the X-fA2 
model did not result in any steric clashes forcing a large 
shift of nonsubstituted atoms during the energy mini- 
mization. The model shows that cGMP bound to the ca- 
nonical ligand-binding pocket, which is highly conserved 
in crystal structures of other cyclic nucleotide-activated 
proteins (Weber and Steitz, 1987; Su et aJ., 1995; Clayton 
et al., 2004; Rehmann et al., 2008), does not tie close 
to the aA-pi region. However, it remains possible that 



%OmM cGMP) 

..A.1..2...3..4...5..6(PB).7..8.B..C 9(3mMcGMP) 

Construct 6 I i.n ^TJ 1 .03 ± 0.05 (7) 

Construct 7 I 
Construct 8 I hi" 



■ " "I 1.52 ±0.18 (9) 

"ill 1.17 ±0.12 (9) 



Construct 9 C 



I 0.26 ± 0.04 (3) 





0.1 1 10 
[Cyclic nucleotide], mM 



R529 

Canonical site 
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the aA-(31 region contributes to a noncanonicaJ cGMP- 
binding site. We emphaisize that although cGMP concen- 
trations >3 mM are required for con action in constructs 
like X-EA.2 with the fCNGA2 BD, this BD also shows gen- 
erally weak binding for cAMP, so that the cAMP dose- 
response curve does not saturate at 3 mM. Thus, the cGMP 
affinity of the binding site responsible for con action 
seems to be comparable to the cAMP afiinity of canonical 
pro-action sites found in the fCNGA2 BD. 

Electrostatic potential calculations on our X-fA2 model 
(Petrey and Honig, 2003) indicate that helix aA and 
strand pi are in, respectively, negative and positive elec- 
trostatic environments (Fig. 4 C) . Although a negative 
environment cannot reliably exclude the binding of phos- 
phate groups (Ledvina et al., 1996) , the positive environ- 
ment around strand pi would seem an intuitive choice 
for binding cGMP For instance, the canonical site has an 
essential arginine R529 (Tibbs et al., 1998), and strand 
pi also contains an arginine R464, which in Construct 8 
was neutralized to glutamine. However, we found that an 
isolated R464Q mutation introduced to X-fA2 did not 
significandy change the low ratio g^so mM cGmp)/§'(3 mM cGMP) 
(Fig. 5 A, X-fA2 R464Q vs. X-£\2) . Similarly, an isolated 
V467T mutation preserved con action and indeed en- 
hanced its observed effects, with a significant decreaise in 
g{30 mM cgmd/S'cs mM cGMP) Compared with that of X-£\2. 
These results argue against strand pi contributing to con 
action at all, and in fact the presence of valine at posi- 
tion 467 in the fCNGA2 BD seems to favor channel open- 
ing compared with the threonine found in rCNGA4. 

In contrast, X-fA2 with an isolated V457E mutation 
shows a mean g(3o mM cGMP)/g(3 mM cGMP) ratio much higher 
than those of the other two point mutants and X-fA2 itself 
(1.26 ± 0.23; w = 7), and significantly greater than unity 
(P < 0.02). The cGMP dose-response curve for X-£\2 
V457E (Fig. 5 B) resembles that of Construct 8 with no 



con action up to 30 mM cGMP. X-fA2 V457E does have 
slighdy higher Pq for a given cGMP concentration than 
does X-fA2, but the increase in efficacy is slight and not 
large enough to raise Pq above 0.5. The low Po, along with 
the significant increase in Pq going from 3 to 30 mM cGMP, 
rule out high Po masking as a cause for the failure to 
observe con action in X-fA2 V457E. This suggests that 
V457 on helix aA is critical for operation of the con- 
action mechanism in CNG channels. 

DISCUSSION 

Our most important new findings in this work are as fol- 
lows. Although there are many BD residues that are not 
conserved between fCNGA2 and rCNGA4, only one is crit- 
ical for con action, namely V457 in helix aA of fCNGA2. 
At least some of the structural determinants of pro and 
con action operate independently of each other: pro- 
action efficacy can be controlled by substituting the PB 
cassette sequence without affecting con-action energet- 
ics, whereas con action can be enabled or disabled by 
substitution of the oA-pi sequence without affecting 
pro action. We produce an example (Construct 3) where 
a con-action mechanism operates efficiendy, yet its sup- 
pressive effects on Pq are partially masked as an indirect 
consequence of efficient pro action. 

Does helix aA contribute to a cGMP-binding site? 

All evidence at hand shows that V457 is critical for con 
action, not simply "contributing" to it. We found no con 
action in any channel that we tested with a V457E muta- 
tion, namely X-rA4, Constructs 6-8, and the point muta- 
tion in X-fA2 background. We emphasize that "critical" 
means "necessary" but does not mean "sufficient"; for 
instance, the presence of valine at position 457 may not 
ensure the operation of con action in every possible 
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sequence context, because we did not find con action in 
Construct 4. We also acknowledge the impossibility of 
confirming the criticality of V457 in all possible sequence 
contexts. Thus, conceivably, con action might be restored 
to X-fA2 V457E by some as yet untested mutations out- 
side position 457. However, any such mutations would 
have to introduce amino acids not found in the fCNGA2 
BD, because X-£\2 V457E currently possesses all fCNGA2 
BD residues except V457. (This contrasts with our suc- 
cessful restoration of con action to Construct 4 by adding 
back noncritical fCNGA2 BD residues to make Construct 
5). Our identification of V457 as critical for con action is 
certainly consistent with known current evidence. Of 
broader importance, the low pro-action efficacy of X-fA2 
V457E (and, more clearly. Construct 8) refutes the suffi- 
ciency hypothesis. Thus, con action does not arise as an 
unavoidable consequence of the rather inefficient pro 
action in fCNGA2; rather, con action can be removed or 
added on top of a preexisting pro-action mechanism by 
virtue of the amino acid at position 457. 

Could V457 then make contact with cGMP, forming 
counterproductive interactions during con action? This 
would constitute a noncanonical cGMP-binding site that 
is structurally distinct from known sites, based on analy- 
sis of our comparative models as follows. We enumerated 
"first neighbor" residues located within 5 A of the resi- 
due 457 side chain in our comparative models (Fig. 6, 
black dashed lines), and then identified "second neigh- 
bor" residues within 5 A of each of the first neighbors' 
side chains (Fig. 6, gray dashed lines) . The residues of the 
"neighbor network" are limited to three noncontiguous 
helices: helix ctA itself, plus helix aD' of the C-linker, 
and helix aB at the C-terminal end of the BD (second 
neighbors only) . All network residues other than V457 
were conserved in all constructs of this study. Based on 
corresponding residues in the HCN2 template structure, 
none of the network residues are proximal to ligand 
bound at the canonical site. Moreover, none of them 
correspond to residues forming the previously reported 
noncanonical cAMP-binding site in Escherichia coli catab- 
olite activator protein (Passner and Steitz, 1997), which 
relies on the ^4-^5 loop along with the DNA BD that 
has no counterpart in CNG channels. 

The neighbor network includes multiple residues 
that could contribute to ligand binding in varied ways. 
Several first neighbors are hydrophobic, suitable for in- 
teraction with the entire valine side chain or with the 
Cp and C-^ methylene groups of glutamate. Additionally, 
however, the salt-bridged residue pair E458+K462 is in- 
cluded, available for stabilizing polar contact with the 
carboxylate group of E457 (H bonds identified by VA- 
DAR software; Willard et al., 2003) . The effect of the 
V457E mutation could thus be a result of any of multiple 
features such as side-chain packing, charge, or specific 
H bonds, and might be as subtle as a slight reorientation 
of one of the three helices. Addressing the large scope 



of mechanistic possibilities will require future mutagen- 
esis in these three helices, which are implicated in con 
action for the first time by our analysis. 

Does helix aA control coupling? 

Another possible role for the Fig. 6 network could be to 
govern the coupling of ligand binding to gate opening, 
like the transmission gearbox in an automobile that 
enables the crankshaft to propel the vehicle in reverse. 
Postulating that the Fig. 6 network is a coupling controller 
and therefore not a noncanonical binding site does not 
force us to assume the existence of a noncanonical site 
somewhere else. Workable con-action mechanisms can 
be conceived in which cGMP binds only to canonical 
binding sites, and it is only the number of vacant versus 
occupied sites that determines whether the coupling re- 
lationship will be reversed (Chan and Young, 2009). 

Helices aD' and ctB in the network have previously 
been identified as coupling-control elements governing 
the efficiency of pro action in CNG channels (Li and 
Lester, 1999; Paoletti et al., 1999; Young and Krougliak, 

2004) ais well as in other cAMP receptors with homologous 
BDs (Rehmann et al., 2003; Zagotta et al., 2003; Kim et al., 

2005) . It may seem odd that these same coupling-control 
elements could reverse the direction of their effects 




Figure 6. Neighbor networli around position 457 proposed by 
comparative modeling. All residues in this structural schematic 
are numbered according to fCNGA2 sequence position. Gluta- 
mate is shown in position 457, with black dashed lines indicating 
its interactions with first neighbors (see Discussion). Parentheses 
indicate residues that are not identified as first neighbors when 
valine occupies position 457. Gray dashed lines indicate interac- 
tions between first neighbors (black type) and second neighbors 
(gray type) . E458 is capable of forming hydrogen bonds or po- 
lar contact with both E457 and K462 (thick dashed lines); other 
neighbor interactions shown are not polar 
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after a certain number of ligands are bound. How- 
ever, occupancy-dependent coupling has already been 
proposed for pro action in CNG channels, based on 
observations that successive ligand-binding events make 
nonidentical energetic contributions even in a homo- 
tetrameric channel (Ruiz and Karpen, 1997; Biskup et al., 
2007) . Ligand-dependent subunit interactions such as 
those mediated by helix aD' (Craven and Zagotta, 2004) 
have been proposed to explain occupancy-dependent 
coupling in pro action (Liu et al., 1998; Ulens and 
Siegelbaum, 2003) as well as the dependence of con ac- 
tion on the arrangement of subunits in a heteromer 
(Chan and Young, 2009) . 

PB cassette region and high Po masking 

Unlike the aA-pi region, the PB cassette contributes 
essential residues to the canonical-binding pocket 
(Altenhofen et al., 1991; Tibbs et al., 1998; Zagotta et al., 
2003; Flynn et al., 2007) , and it is therefore not surprising 
that significant effects on pro-action efficacy result from 
substitution of the fCNGA2 caissette. We caution that the 
effects of PB cassette substitution are expected to depend 
on the background sequence context, as in any structure- 
flinction study. Unfortunately, the structural effects of PB 
cassette substitution cannot be reliably explored through 
our comparative model, because fCNGA2 contains a 
six-residue insertion in its cassette region with no ho- 
mologous counterpart in the cassette of HCN2. Further 
investigations of efficacy determinants may benefit from 
testing cassettes from different CNG channels, because 
the insertion is shorter in some CNG family members. 

The broader significance of the PB cassette replace- 
ment producing Construct 3 is that it illustrates the 
principle of masking through an increased gating equi- 
librium constant. The persistence of efficient con ac- 
tion in Construct 3 is detectable when this masking is 
taken into account, refuting the necessity hypothesis. 
Because low pro-action efficacy is not a prerequisite for 
the operation of con action, it may be that the con-action 
mechanism exists in some other previously studied CNG 
channels but has gone undetected because of masking 
(see discussion of rat CNGA2 in Biskup et al., 2007, and 
Chan and Young, 2009). 
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